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Abstract A trapped Bose-Einstein condensate, being strongly perturbed, exhibits several 
spatial structures. First, there appear quantum vortices. Increasing the amount of the injected 
energy leads to the formation of vortex tangles representing quantum vortex turbulence. 
Continuing energy injection makes the system so strongly perturbed that vortices become 
destroyed and there develops another kind of spatial structures with essentially heteroge¬ 
neous spatial density. These structures consist of high-density droplets, or grains, surrounded 
by the regions of low density. The droplets are randomly distributed in space, where they can 
move; however they live sufficiently long time to be treated as a type of metastable creatures. 
Such structures have been observed in nonequilibrium trapped Bose gases of *^Rb subject 
to the action of an oscillatory perturbation modulating the trapping potential. Perturbing 
the system even stronger transforms the droplet structure into wave turbulence, where Bose 
condensate is destroyed. Numerical simulations are in good agreement with experimental 
observations. 

PACS numbers: 03.75.-b, 03.75.Hh, 05.30.Jp, 47.27.E-, 47.27.Gs, 67.40.Vs 


1 Introduction 

Bose atoms in a trap can be cooled down to such low temperatures when almost all of them 
pile down to the ground state corresponding to Bose-Einstein condensate. If the trapped 
Bose-condensed cloud is slightly perturbed, it exhibits elementary collective excitations, 
similarly to other finite systems, such as quantum dots, atomic nuclei, or clusters mi2ll3ll4l 
[^. But what happens when the Bose-condensed cloud is strongly perturbed? 
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We know that, when the external perturbation imposes a rotation moment, there can 
arise quantum vortices that may form a vortex lattice EB. But if the external perturbation 
imposes no rotation moment, but just strongly shakes the condensate, as has been suggested 
in Refs, mi, then no vortex lattices arise. Instead, the increasing number of vortices forms a 
random tangle corresponding to quantum turbulence lll0lllllll2lfT3l . This vortex turbulence 
of trapped atoms has been observed in three-dimensional traps ll4IIL5lll6ll as well as in 
quasi-two-dimensional traps (13, and is summarized in the review articles (unsiiaoi. 

In the present paper, we study what could be other possible nonequilibrium states in a 
trapped strongly perturbed Bose-Einstein condensate. We investigate this problem both in 
experiments with trapped atoms of *^Rb and accomplishing computer simulations for the 
setup exactly corresponding to these experiments. Observed experimental data are in good 
agreement with numerical simulations. 


2 Methods of strong perturbation 

Generally, a Bose-Einstein condensate is characterized by an equation for the condensate 
wave function rj =rj(r,t). For a system of weakly interacting atoms at temperatures close to 
zero, the condensate-function equation takes the form of the nonlinear Schrddinger equation 
(NFS) 

Here and in what follows, the Planck constant is set to one, U = U (r,r) is the total external 
potential, and the atomic interaction strength is 

0 = 471 — , ( 2 ) 

m 

with flj being scattering length and m, atomic mass. 

When the system is in equilibrium, the condensate is in its ground state. To perturb the 
condensate, it is necessary to inject energy into the system. The amount of injected energy, 
per atom, during the excitation time f, can be represented as 



where N is the total number of atoms, and the energy Hamiltonian is 

H = I ri{r,t) (^-^ + u'^ri{r,t)dr+^^ j \ r\{r,t) \‘^ dr . (4) 

This tells us that energy can be injected through one of two ways, either adding to 
the trapping potential U{r) a time-dependent perturbative term, so that the total external 
potential be the sum 

U{r,t) = U{r)+V{r,t), (5) 

or by varying the scattering length, so that the interaction strength be time-dependent, 

0{t)=4n‘^^. ( 6 ) 

m 

The perturbations are assumed to be organized in such a way that no rotation moment is 
imposed onto the system. 
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If the trap modulation is done by means of an alternating potential 


V(r,t) ~Acos(®t) (A > 0) , 


as is advanced in Refs. mil, then the injected energy, for sufficiently long time t ^ iTtjw, 
takes the form 


Ein i ~ — A(Ot , 

K 

which makes it possible to study the amplitude-time phase diagram through the relation 


(7) 


A 


TZEifi j 

2cot 


( 8 ) 


The other possibility of injecting energy is by modulating the scattering length, making 
it time-dependent with the help of Feshbach resonance techniques, as has been suggested 
in Refs. I21II22II23I . The time dependence of the scattering length is realized through the 
oscillating magnetic field, when 


as{t) =as 


1 - 


AB 




(9) 


where is a background scattering length, AB, resonance width, and B„ is a resonance 
field. The modulating magnetic field oscillates by the law 

B{t) = Bq +Bi cos(a)t) . (10) 

For a small oscillation amplitude, such that |5i/5o| ^ 1, the effective scattering length is 


as{t) c^OQ+ai cos(a)t) 


with the notation 


ao = aA 1 - 


AB 


Cl\ = Gs 


BiAB 


Bo-B„ 

Then the interaction strength takes the form 

(p{t) = 00 + ^’i cos(a)r) 

in which 


(So-5„)2 ■ 


00 = 4;i' 


ao 


01 = 4;r 


Gl 


The injected energy reads as 


^inj ~ ^ p 01 cot , 


where p is average atomic density. Denoting the effective amplitude 


( 11 ) 


( 12 ) 


(13) 


1 fli 

A = - p0i = 2;rp — , 

2 m 

we get the same amplitude-time relation (8) as in the case of the trap modulation. Thus, these 
two methods of perturbing condensate should lead to the same consequences, producing 
similar strongly excited nonequilibrium states. 
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The nature of the possible nonequilibrium states essentially depends on the atomic scat¬ 
tering length as well as on trap geometry. Important quantities quantifying these character¬ 
istics are the aspect ratio 



(14) 


in which I±_ 
parameter 


1 jlz = ^ly/fncOz are the oscillator lengths, and the effective coupling 


a, 

g = 47tN — . 


(15) 


It is possible to use other methods for perturbing Bose-Einstein condensate, for instance, 
by creating laser-induced obstacles (24l, which is equivalent to the perturbation by external 
potentials, that is, to trap modulation. 


3 Nonequilibrium atomic states 

The system is initially prepared in an equilibrium state, at temperature close to zero, with 
almost all atoms being in Bose-Einstein condensate. Energy is injected into the system by 
modulating the trapping potential as described in Refs. fnn^uToinm . Injecting energy into 
the trap produces a sequence of nonequilibrium states that we have studied experimentally 
as well as by computer simulations for the order-function equation (1), adding there the 
attenuation by replacing in Eq. (1) i by i — y, which imitates the loss of energy and atoms 
(una. The experimental results are in very good agreement with computer simulations. 
Physical effects, occurring during the process of perturbing the system can be illustrated by 
theoretical estimates that are in close agreement with both experimental results and numeri¬ 
cal simulations. 

Weak perturbation. Starting pumping into the trap the injected energy (3), first, creates 
elementary collective excitations corresponding to small fluctuations around the equilibrium 
state. This regime of weak perturbation lasts in the interval of energies 

0 < Ei„j < 2E,,or , (16) 

where the injected energy is smaller than the energy required for exciting at least a pair 
of vortices with opposite vorticities. Such a pair makes the total vorticity of the atomic 
cloud zero, as it should be, when the perturbing potential does not impose rotation. It is 
worth stressing that the created pair does not form a bound state, but the vortices are rather 
independent of each other. The energy of a vortex per atom can be defined I25II26I as 

= (17) 

Keeping in mind the setup of experiments I14II15II16II201 with ^^Rb, we have the trap 
characteristics 

cox = 1-32 X 10^s^', 0)^ = 1.45 X lO^s^' , 

= 0.74 X 10^"^cm , 4 = 2.25 x lO^'^cm {^'^Rb) , 

which implies the aspect ratio a = 0.11. The *^Rb atoms have the scattering length a^ = 
0.577 X 10^® cm. With the number of atoms N Kilx 10^, the effective coupling (15) is g = 
1.96 X 10^. This gives the pair vortex energy as lE^or = 0.566 x 10^*^ eV. It is convenient 
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to measure the energies in units of the transverse trap energy Ex_ = h(j)± that in the present 
case is E±_ = 0.869 x 10^*^ eV. Then 2Ey„r = 0.65Ztx- 

In order to stress that the appearance of vortices as well as of other nonequilibrium states 
essentially depends on the involved sort of atoms and on the trap geometry, we compare the 
typical values for atoms with the setup employed in the group of Hulet I27II28I dealing 
with ^Li atoms and a more elongated trap, where 

= 1.48x lO^s-', = 0.304 X lO^s-', 

= 2.5 X 10^"^cm , Zj = 1.7 X lO^^cm CLi) , 

which gives the aspect ratio a = 0.021, an order smaller than in the case of ^^Rb. The 
scattering length of ^Li is tuned by Feshbach resonance to = 3.2 x 10^^ cm, which is 
two orders shorter than that of Rb. The effective coupling g = 0.48 x 10^ is much smaller 
than for ^^Rb. The energy, required for vortex generation, is 2E^,or = 1.45 x 10^^^ eV. The 
transverse energy here is E± = 0.97 x 10^*^ eV. Hence 2£,,or = 0.49E±. In this setup, it 
is more difficult to create vortices, since the required injected energy Ei„j should be about 
three times larger than in the case of ^^Rb. 

In the regime of weak perturbation (16), vortices and other topological modes are not 
created, unless the modulating potential is tuned to a transition frequency a)„ = E„ — Eq, 
corresponding to a transition between the ground-state condensate, with energy Eq, and a 
coherent topological mode t8l l20ll22ll23ll291 . with energy £„. The condition ® = ®„, when 
topological modes are generated, even under weak perturbation, can be called topological 
resonance. Note that this is not the standard parametric resonance (^, when O) = 2cOq, 
where coq is a system natural frequency, or more generally, when kco = 2coo, where k = 
1,2,.... The natural frequency coq characterizes the motion of the system as a whole and, for 
trapped atoms, is defined by the effective oscillator trap frequency. However the transition 
frequencies co„ describe the internal transitions between the collective energy levels and, for 
an interacting system, (0„ can be very different from the effective oscillator trap frequency 
[ 3 TII 32 I . Topological modes can also be created under harmonic generation, when kco = C0„, 
where k = 1,2,..., under parametric conversion, when there are two external modulating 
fields, with the frequencies cOi and CO 2 such that cOi ± CO 2 = C0„, and under combinations of 
these resonance conditions. But if there is no a kind of a resonant tuning, the topological 
modes do not arise. 

When considering nonequilibrium states, it is necessary to distinguish the atomic cloud 
itself and the subsystem of its elementary excitations. The atomic cloud can be slightly per¬ 
turbed, but, nevertheless, such a perturbation can result in the creation of many elementary 
excitations. The waves, corresponding to these excitations, are involved in wave interac¬ 
tions and are described by quantum kinetic equations 11331134113511 . In the subsystem of these 
excitations, there can arise weak wave turbulence I36II37L which is inhomogeneous in the 
presence of a trap 1381 . 

Vortex formation. As soon as the injected energy surpasses the vortex energy (17), 
vortices start being created. The regime of vortex formation corresponds to the interval 

2E,,or < Einj < E,ur , (18) 

until the number of vortices is so high that the system turns into turbulent regime, which 
happens when the number of vortices reaches a critical value Nc, corresponding to the energy 

Tt'tur — blcTSyor ■ 


(19) 
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The critical number can be defined by the condition that the mean distance between the outer 
parts of the vortices becomes equal to their diameter 2^, given by twice the coherence length 
~ fijmc, where c ~ (h/m)^^4Kpas is the sound velocity. This means that the distance 
between the centers of the vortices is 4i§, which gives 

( 20 ) 

For the experiments with ^^Rb, the atomic cloud radius and length are rx w 4 x 10^^ cm and 
L w 6 X 10^^ cm, respectively. The coherence length is w 2 x 10^^ cm. Thus, the critical 
vortex number is A'c ~ 25, which is in perfect agreement with experiments and simulations. 
This corresponds to the vortex turbulence energy Emr = 0.71 x 10^^* eV, or in units of the 
trap energy, E,ur = ?>.2E^. 

It is much more difficult, if possible at al., to create vortex turbulence in an elongated 
trap with trapped ^Li, as in the experiments of the Hulet group I27II28I . In the latter case, 
the cloud radius and length are ~ 3 x 10^^ cm and L fa 2 x 10^^ cm, while the coherence 
length is ft! 2 X 10^^ cm, so that a single vortex would fill almost the whole trap. Formally, 
the critical number (20) would be 0.14. However, there is no turbulence with a single vortex. 

Typical vortex states are illustrated in Fig. 1 showing the results of numerical calcula¬ 
tions for the atomic density in a transverse cross-section of the trap, for the setup character¬ 
ized by the experimental data with ^^Rb atoms. The higher density corresponds to brighter 
colour. So that vortices are shown as black spots. By calculating the vorticity, we find that 
practically all vortices have vorticity plus or minus one. Typical vortex states observed in 
the time-of-flight experiment are presented in Fig. 2. Again, brighter colour corresponds to 
higher density. 



Fig. 1 (Color online) Vortex state. Typical distributions of the density of atoms in a transverse cross- 
section of the trap, found by numerical simulations. Brighter colour corresponds to higher density. So that 
vortices are seen as black spots. 


Increasing the amount of the injected energy produces a larger number of vortices. Since, 
as is clear from relation (7), the injected energy is proportional to the modulation time, the 
number of vortices, in the vortex state, increases with time. After their number reaches 
the critical number 25, the regime of vortex turbulence comes into play, when the number 
of vortices yet increases, but slower than before. This behavior is illustrated by numerical 
simulations presented in Fig. 3, showing the number of vortices as a function of time, under 
fixed amplitude. The results are in good agreement with experimental data. We see that at 
some moment of time the number of vortices abruptly diminishes. Then the system passes 
to another regime, termed grain turbulence or droplet turbulence, to be described below. 
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Fig. 2 (Color online) Vortex state. Experimental observation of vortices in a transverse cross-section of the 
cloud of atoms in time-of-flight measurements. Brighter colour corresponds to higher density. So that 
vortices are seen as black spots. 



Fig. 3 (Color online) Number of vortices, found by numerical simulations, as a function of time, that is, as a 
function of injected energy. 


Vortex turbulence. The regime of vortex turbulence lasts till the number of vortices 
becomes sufficiently high, so that the mean distance between vortices is such that their 
interactions become rather strong, when the interaction energy of two vortices at the distance 
5 from each other becomes equal to the vortex energy. Equating these energies, given in Ref. 
(26l, we get 


In 




where R is an effective radius of the cloud. This yields 


5 = 1.19 


(21) 


Taking the Thomas-Fermi radius 



R = l 


1/5 
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results in the critical distance 


5 = 0.86 {agf 


10 


( 22 ) 


For ®^Rb, the latter is 5 ~ 0.71 x 10^“* cm, so that 5 ~ 3.6 ^. 

One distinguishes random vortex turbulence, or Vinen turbulence, as opposed to the 
correlated vortex turbulence, or Kolmogorov turbulence I19II39II40II41II42I . In Vinen tan¬ 
gles, far-field effects tend to cancel out, the motion of a vortex line is mainly determined by 
its local curvature, and the vortex tangle is homogeneous. The intervortex distance in Vinen 
turbulence is much larger than the coherence length. Therefore the summary energy of vor¬ 
tices in the random Vinen tangle can be approximated by the product of the vortex number 
and a single vortex energy. 

In the regime of vortex turbulence, we have observed, the distance between the vortices 
is yet such that the energy of vortex interactions is smaller than the energy of a vortex, 
although interactions can be noticeable. The vortices are not correlated by the used trap 
modulation that does not impose such correlations, nor imposes a rotational polarization. 
However, to distinguish precisely what type of turbulence has been realized requires the 
study of energy spectra, which is yet in progress. In the present publication, we concentrate 
on the study of density distributions that can be directly observed and which allows for the 
qualitative classification of the observed nonequilibrium regimes. 

The regime of vortex turbulence is realized in the interval of the injected energy 


Etur < Ei„j < Efog , (23) 

where the upper energy boundary 

Ef„g=N:E,or (24) 

corresponds to the critical number of vortices, when interactions become rather strong, 
which means that the distance between the vortices is smaller than the critical distance (22). 
This critical number can be estimated as 



For the considered case of ^^Rb, we have N* = \.2Nc w 30. This defines the energy Efog = 
0.84 X 10^** eV, corresponding to ^fog — ‘^eie^. 

It is clear that the estimates for the energies Etur and Efog give their approximate values. 
It turns out, these values are close to each other. Although the absolute range between the 
energies Etur and Efog, of course, is approximate, but its narrowness means that the region 
of vortex turbulence is quite narrow. This has really been observed in experiments as well 
as in numerical modelling. 

The typical cross-section of the random vortex tangle, obtained in numerical simula¬ 
tions, is demonstrated in Fig. 4. And Fig. 5 shows the vortex turbulent state observed in the 
experiment. 

Grain turbulence. After the number of vortices in the random tangle grows to the crit¬ 
ical number N*, the vortices start strongly interacting with each other by colliding and de¬ 
stroying each other. Their number sharply decreases, as is shown in Fig. 3. The remnants 
of the destroyed vortices form the pieces of the granulated Bose condensate, which can be 
called droplets or grains. These droplets, being multiscale in the interval (1 — 5) x 10^^ cm, 
have the typical size ~ 3 x 10^^ cm, which is close to the coherence length ^ ~ as it 
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Fig. 4 (Color online) Vortex turbulence. Cross-section of a random vortex tangle obtained by numerical 
simulations. Darker coulor corresponds to lower density, thus vortices are shown as black spots. 



Fig. 5 (Color online) Vortex turbulence. Experimental obseiwation of a vortex turbulent state in the cloud 
*^Rb atoms. Multiple dai'k spots correspond to vortices. 


should be for a coherent object. The coherence of a droplet is confirmed by calculating its 
phase that is constant through the droplet volume. The droplets are surrounded by a rarefied 
gas of much lower density; the ratio of the droplet density to the density of surrounding 
Pdr/Psur Varies in a wide range: from a few tens to 100. At each snapshot, the droplets are 
randomly distributed in space, forming no spatial structures, reminding water droplets in 
fog. Therefore, this regime can be termed droplet turbulence or grain turbulence l43l . The 
droplets are compact objects, with the linear sizes in different directions being close to each 
other. The typical droplet lifetime is of order 10^^ s. The droplets do not form any regular 
structure, since the average atomic density is not high, being much lower than that necessary 
for creating a kind of crystalline order 144II45I . The random spatio-temporal distribution of 
the droplets makes it admissible to call this state as grain turbulence. 

In numerical simulations, considering the opposite process of Bose condensate forma¬ 
tion from an uncondensed gas, the regime of grain turbulence corresponds to that of strong 
turbulence II46II471 . We prefer to call it grain turbulence, since this term reflects the physical 
nature of the system consisting of dense droplets, or grains, in a rarefied surrounding. 

The equivalent interpretation of the granular state could be by imagining it as a col¬ 
lection of dark solitons, typical of defocusing NLS, inside a Bose-condensed surrounding. 
The granular state is similar to the critical balanced state that is a state where the linear and 
nonlinear terms are balanced for a wide range of scales mm- 

It is worth emphasizing that the transition from one dynamical regime to another is not 
an abrupt phase transition, but a gradual crossover, although it can be very sharp. Therefore 
the transition lines are, of course, conditional, showing where one dynamic behavior changes 
to the other. At the same time, some features of one dynamic regime can survive to the 
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Fig. 6 (Color online) Grain turbulence. Cross-section of the atomic cloud composed of dense Bose- 
condensed grains, randomly distributed in space, inside a rarefied surrounding. Droplets are seen as bright 
spots. 



Fig. 7 (Color online) Grain turbulence. Experimental observation of a granular turbulent state in the cloud of 
*^Rb atoms. Brighter spots correspond to higher density, while darker regions, to rarefied density. Droplets 
are seen as bright spots inside a rarefied surrounding in dark. 


Other region. For instance, in the regime of grain turbulence, some occasional vortices can 
occur, although their number is very small and does not change the overall picture. However, 
eventually a few vortices can arise and then annihilate in various locations at random. 

Continuing pumping energy into the trapped system destroys Bose condensed droplets. 
So that the grain turbulence exists in the interval of energies 

Efog<Einj <Ec , (26) 

until the injected energy per atom is so high that all Bose condensate becomes destroyed, 
when the energy reaches the value 

Ec = kBTc, (27) 

where Tc is the critical temperature of Bose-Einstein condensation. For the treated case of 
®^Rb, the critical temperature is Tc = 2.76 x 10^^ K. Then we have Ec = 0.238 x 10^'° eV. 
This gives Ec = 21AE^. 

Figure 6 presents the cross-section of the atomic cloud in the regime of grain turbulence, 
found numerically. Droplets are shown as bright spots surrounded by dark rarefied regions. 
Each droplet is formed by Bose-condensed atoms, which is confirmed by a constant phase 
inside a droplet. In Fig. 7, an experimentally observed granular state is shown. 

Wave turbulence. Increasing the amount of injected energy should finally lead to a 
complete destruction of condensate (SOj. This happens after the injected energy surpasses 
the critical energy Ec, 

Ein J A Ec , 


(28) 
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Fig. 8 (Color online) Wave turbulence. Cross-section of the trapped cloud in the regime of wave turbulence, 
found numerically. The left figure shows the density distribution, where brighter colour corresponds to higher 
density. The right figure shows the spatial phase distribution, where the larger phase is of brighter colour. 


Grain 

turbulence 

Wave 

turbulence 

PdJ Ps.r ~ 100 

P^IPsur ~ 3 

constant phase inside 

random phase inside 

^kin ! ^int ~ ^ 

^kin ! ^int ~ 300 

~ 10^^ cm 

~ lO^^'cm 

hn - 1 

L » ^ 


Fig. 9 (Color online) Characteristic features of the regimes of grain turbulence and wave turbulence. 


when all Bose-Einstein condensate is getting destroyed. Then the system is represented by 
small-amplitude waves, because of which it is termed wave turbulence, or weak turbulence. 
The linear sizes of waves are in the interval (0.5 — 1.5) x 10^“* cm, so that the typical wave 
size is l„ ~ 10^^ cm. The wave density is only slightly greater than that of their surrounding, 
Pw/Psur ~ 3. The phase inside a wave, as well as between the waves, is completely random. 
The system kinetic energy is much higher than the interaction energy, so that the system 
can be represented as a collection of almost independent modes described by a kinetic equa¬ 
tion with four-wave processes 02). The dynamic transition from grain turbulence to wave 
turbulence is a crossover, such that in a wide region grains coexist with waves. 

Also, in the weak turbulence regime, there exist the so-called ghost vortices that are 
nodal points of the wave field, which are present in abundance if the waves are linear OH- 
The cross-section of the trapped cloud in the regime of wave turbulence is shown in 
Fig. 8, where the system density and phase are presented. This figure is obtained by numer¬ 
ical simulations. Because of the large amount of injected energy, required for creating this 
regime, it has not yet been reached in experiments. 

To stress the principal difference between the regimes of grain turbulence and wave 
turbulence, their characteristic features are compared in Fig. 9. 

In conclusion, we have studied both experimentally and by means of computer simu¬ 
lations strongly nonequilibrium states of trapped atoms, prepared in Bose-condensed 
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Fig. 10 (Color online) Sequence of nonequilibrium states produced by modulating the trapping potential for 
*^Rb atoms. 



Fig. 11 (Color online) Numerical amplitude-time phase diagram for *^Rb. 


State and subject to strong perturbations realized by modulating the trapping potential. The 
sequence of the nonequilibrium regimes, obtained by computer simulations, is summarized 
in Fig. 10. And the numerical amplitude-time diagram is presented in Fig. 11. The latter is 
in good agreement with the experimentally found amplitude-time diagram discussed earlier 
I15II16I . It is interesting to note that the sequence of the nonequilibrium states, generated 
in the process of strong perturbation, which includes weakly nonequilibrium state, vortex 
state, vortex turbulence, grain turbulence, and wave turbulence, repeats the analogous states, 
although in the reverse order, as those occurring during the process of equilibration I51II52I 
through a nonequilibrium Bose-condensation phase transition from a strongly nonequilib¬ 
rium uncondensed state to an equilibrium condensed state I37II46II47II53I . Thus, we generate 
the states typical of the Kibble-Zurek mechanism I54II55I . but in the opposite temporal or- 
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der. More discussions on such an inverse Kibble-Zurek scenario will be given in a separate 
paper. 
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